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A self-consistent, spin rotational invariant Green's function procedure has been devel- 
oped to calculate the spectral function of carrier excitations in the spin-fermion model 
for the Cu02 plane. We start from the mean field description of a spin polaron in the 
Mori-Zwanzig projection method. In order to determine the spin polaron lifetime in the 
self-consistent Born approximation, the self-energy is expressed by an irreducible Green's 
function. Both, spin polaron and bare hole spectral functions are calculated. The numer- 
ical results show a well pronounced quasiparticle peak near the bottom of the dispersion 
at (7r/2, 7r/2), the absence of the quasiparticle at the F-point, a rather large damping 
away from the minimum and an asymmetry of the spectral function with respect to the 
antiferromagnetic Brillouin zone. These findings are in qualitative agreement with pho- 
toemission data for undoped cuprates. The direct oxygen-oxygen hopping is responsible 
for a more isotropic minimum at (ti /2,n /2). 

71.27.+a, 74.25.Jb, 74.72.-h 



1. INTRODUCTION 



In recent years there was some progress to answer the important question on the nature of quasiparticlcs 
in cuprate superconductors, but the task is not yet completed. Usually, the problem is reduced to that 
of thC|-halc motion in a two dimensional quantum antifcrromagnet. Many work was devoted to the t-J 
modelETtl treating the antiferromagnet as a two-sublattice Neel state. It was shown that the hole motipn 
is characterized by a quasiparticle band of bandwidth ~ 2 J that is separated from the incoherent part.EFEl 
The minimum of the dispersion is found at the momentum (7r/2,7r/2) (lattice constant g — 1) nearly 
degenerate with the energy at (ir, 0). Such a result was obtained by several methods. One of them is the 
self-consistent Born approximation (SCBA) for a reduced version of the t-J model consisting of spinless 
holes coupled to spmons.B'ffirLl That method became especially popular since it is sufficiently simple and 
allows the calculation of the spectral function. It gives rise to a dispersion and also a spectral function 
which is symmetric jsrith respect to the magnetic Brillouin zone (BZ). Similar results were obtained for 
the Emery model. Bil3 _ . 

In a remarkable experimentO the quasiparticle dispersion of such a spin polaron was measured di- 
rectly by creating a photohole in the antiferromagnetic insulator Sr^CuC^C^. Whereas the width of the 
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dispersion was in reasonable agreement with the theoretical predictions, its shape had some deviations 
especially at the point (n, 0) in momentum space. Afterwards it was shown that the deviations can. be 
diminished by taking into account hopping terms to second and third neighbors in the t-J model.E3ii.il 
These additional terms correspond to a proper treatment of direct oxygen-oxygen hopping in the Emery 
model.E3o However, some deviations between theory and experiment remain. So, the measurement 
shows no quasiparticle weight at the T-point (0, 0) but the spin polaron in the t-J model has. Next, 
one observes quite a remarkable damping especially away from the minimum. Some theories (especially 
the usual SCBA) are symmetric with respect to the antiferromagnetic (AFM) BZ, but the experimental 
data show a strong asymmetry: after crossing the AFM BZ boundary (coming from the T-point) the 
quasiparticle looses its weight nearly abruptly. That .deficiency, however, seems to be merely a property 



of the usual SCBA and not of the t-J model itsclfJljl 

In the usual SCBA there is no dispersion of the bare spinless hole; the dispersion appears only due to 
the coupling to the spinons. Therefore, one might expect large vertex corrections. In fact, however, the 
low order vertex corrections cancel in the pure t-J model as it was noted in Refs. ||,^| Such a cancelation 
is absent in the more realistic Emery model. But there are many indications (see also Ref. [l7]), and we 
will add some arguments below, that the photoemission data cannot be explained by the t-J model alone, 
at least not throughout the whole BZ. 

We concentrate-in the present paper on the spin-fermion model (the reduced form of the three-band 
Emery model)ll3'E3 for the Cu02 planej-apd present a new scheme to calculate the spectral function 
starting from the mean field descriptionfSEil of the spin polaron. One can expect that the fluctuations 
around the mean field spin polaron are smaller than around the dispersionless spinless hole. The mean 
field description of the spin polaron in the Emery model was introduced in Refs. |2^,^l],^6|. It is based on 
a description of the hole motion not in the two-sublattice Neel state but in the spherically symmetric (in 
spin space) AFM state with long range order (LRO). Such a spin singlet state is characterized here by a 
vanishing expectation value of the magnetization (Sf ) = (the only possible nonzero expectation values 
are scalars with respect to rotations in spin space), but a finite value of the spin-spin correlations at 

infinite distance (SoSr) — ► (— 1) R M for R — > oo. It can be constructed, for instance, from a resonating 

valence bond (RVB) state with long-range bonds that decay not too fast.^l The spherically symmetric 
state is an appealing candidate to study the transition from the paramagnetic spin liquid (with only short 
range correlations) to a state with long-range order. 

A spherically symmetric Green's function (GF) theory for Heisenberg models was constructed by Kondo 
and Yamaji (KY)E3 and it will be used as one ingredient of the present theory. It is a GF decoupling 
scheme which is applicable to cases in which there is no symmetry breaking, especially to low dimensional 
systems (for the 2D square lattice see Refs. |2^^). As a result, the carrier excitations (moving in the 
spherically symmetric background) preserve the periodicity relative to the full BZ (not the magnetic one). 
Namely this leads to a qualitative agreement with the experimental results in the undoped compound 
Sr2Cu02 CI2. An explanation might be that the actual experiment was done at room temperature, 
whereas the Neel temperature of Sr2Cu02Ci2 is only 250 K. 

The aim of the present work is twofold. In the first place it is a methodical work proposing a new 
technique to calculate the spectral function of quasiparticles with complicated operator structure. Our 
technique combines the Mori-Zwanzig projection method with a self-consistent calculation of the .self-, 
energy in terms of irreducible GF (Tserkovnikov techniqueEH) . In difference to our recent worksE3E3 
we introduce a simultaneous description of a spin polaron and a bare hole, and investigate the effects 
of direct oxygen-oxygen hopping which only gives us the possibility to reach a realistic description. 
Furthermore, we discuss also the upper parts of the spectrum in the spin-fermion model whose lowest 
branch corresponds to the spin polaron. The second aim is to present numerical results for the spectral 
function. We will distinguish three regions in k-space where the spin polaron exists as a quasiparticle 
with infinite lifetime, with finite lifetime or where it is strongly overdamped. Our method leads to several, 
qualitatively new features which are not present in the SCBA of spinless holes coupled to spinons. But 
these features may be observed in the experimental data for Sr2Cu02Cl2- We are concentrating here on 
a qualitative discussion rather than fine tuning the parameters of the Emery model in its spin-fermion 
form to obtain maximal agreement with experiment. 

Our paper is organized as follows: after presenting the set of basis operators which define the mean 
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field description of the spin polaron (Sec. 3) we show the general procedure to calculate the damping 
(Sec. 4) and calculate the corresponding vertex (Sec. 5). The self-consistent set of equations for the GF 
is given in Sec. 6 and the results are presented in Sec. 7. 



2. SPIN-FERMION MODEL 

An extra hole propagating in the Cu02 plane of cuprates strongly interacts with the antiferromag- 
neticalhixorrelated copper spin subsystem. The main features of the hole motion are described by the 
modebm 

H = T + J + h, (2.1) 
r = ^rJ2pl(l+Sn)m, J = ^Y, S * S *+ S > ( 2 - 2 ) 



2 "i"v , 2 
R R,g 



■£[4, 

R 



(cR+a y + CR-a s + C R+gx+Bly + C R+gj _ a J + h.C. , (2.3) 



with 

PR 



g 



where a x>y = ^g x . v ,S = ±g x ,±Ey 

Here and below the summation over repeated indexes is understood everywhere; {...,...}, [...,...] 
stand for anticommutator and commutator respectively; g x . a are basis vectors of a copper square lattice 

( |g| = 1), R + a are four vectors of O sites nearest to the Cu site R; the operator c R+a creates a hole 
predominantly at the O site (the spin index is dropped in order to simplify the notations); <r Q are the 
Pauli matrices; the operator S represents the localized spin on the copper site. We do not introduce the 
explicit relative phases of p- and d -orbitals since they can be transformed out by redefining the operators 
with phase factors exp (iqoR) . qo = (tt, tt). In order to compare our results with other authors and with 
experiment we restore these phases by changing in the final results k — > k' = k — q . The parameter r is 
the hopping amplitude of oxygen holes that take into account the coupling of the hole motion with the 
copper spin subsystem, J is the constant of nearest neighbor AFM exchange between the copper spins. 
For the parameter values we take r ~ 0.5eV, J = 0.2t throughout the paper, h — 0, or h = 0.3t to 
illustrate the influence of direct oxygen-oxygen hopping. 



The Hamiltonian (2.1) corresponds to the regime Ud ^S> \s p — Ed\ ^> t v d in the Emery model. It 



can be obtainc dEl e. g. by a canonical transformation of bare fermion operators a — e s ae s . Then the 



transformed operators a which enter the Hamiltonian (2.1) acquire the admixture of bare operators at 



adjacent sites. For example, the "oxygen hole" operator in the Eq. (2J3) may be expressed as 



CR+a.s = C S C R+a!S e 5 W C R+SL , S + -. — rG?R. s ( 1 - d|l - s rf R.-s ) ■ (2-5) 

This operator annihilates the hole in a state that is an antibonding (in hole notation) combination of the 
oxygen p„ state and the singly occupied copper d x 2_ y 2 state at adjacent sites. 

The non-bosonic character of spin operators and the com plexity of the spin subsystem ground state 



| Q) make a perturbative treatment of the Hamiltonian (2.1) difficult. For this reason we develop in the 
next sections pthe-imethod of Green's function calculation which combine the Mori-Zwanzig projection 
method (PM)E!rE3 with the Tserkovnikov technique. The PM provides the possibility to allow for the 
local correlations thoroughly. For example, the local constraints are naturally included in tin-calculation. 
We thus obtain the proper description of small polaron (so called Zhang-Rice (ZR) singlet formation 
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that is governed by the Kondo term in f . On the other hand, the PM always provides the Green's function 
as the sum of simple poles, since only a finite number of basis operators may be included in the calculation. 
The consideration of polaron scattering on spin excitations, which is responsible for the damping of the 
polaron state, demands an infinite extension of the basis set. This can be done effectively by the self- 
consistent -Born approximation (SCBA) for the polaron self-energy, calculated within the Tscrkovnikov 
technique.c3 That means that we suppose the effect of the residual polaron-spin interaction to be rather 
small and tractable within a self-consistent perturbative scheme. The main point here is the appropriate 
choice of the operator basis set for PM and a careful division of the interaction into mean field and 
residual part. 



3. BASIS OPERATORS AND PROJECTION METHOD 



^From the very beginning we want to take into account properly the local correlations f without the 
violation of spin commutation relations. We consider a spin-liquid state with spin rotational symmetry 
for the copper subsystem, e.g. the spin-spin correlation functions satisfy the relation C r = (SrSr+j.) = 



3 Si 



x{y,z) gx(y,z) 



R 



R+r 



We choose a set of three basis operators. The first two of them constitute the ZR 



polaron 



1 



R 



-ikR 



PR , B 



1 



2,k 



Vk>/N 



J>- lkR 5RPR, 



(3.1) 



R 



where the factors /3k and arise due to the orthonormalization 



/?k = V 1 + 7k , fk = 



Cg7k 



(3.2) 



and 7k = \ (cos k x + cos k y ). In fact, the basis operators ([O]) can be combined in two ways corresponding 
to the ZR singlet and part of the ZR triplet state. Only the lower lying singlet combination builds the 
spin polaron quasiparticle. We can also write 



where 



ly\ b — — 

' /?k 



cos [ — I c k , 



COS \ — ) Ck,„ 



(3.3) 



1 -«k(R+a x „) 

7W ? e R+s 



Now it is easy to see that one has to introduce the operator 

' k, 



B. 



3,k 



1 



Ck,; 



Ck,-, 



(3.4) 



in order to represent the full set of bare hole operators. That is important if we want to find the hole 
spectral weight of the spin polaron band. 

(i) 

In a first step we construct the three eigenoperators Bi.k and the corresponding bands fly. in the mean 
field approach: 



Bi,k = af (k) B j>k 



£ k = Bi, k 



(i) 



(3.5) 



where we introduced a simplified notation for the most interesting lowest spin polaron band _Qi_with 
the polaron annihilation eigenoperator B^. We use the projection method of Mori and ZwanzigpllEl and 
introduce the retarded Green's functions (GF): 
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^ dte^-t'HiB^it)^]^')}) . 



(3.6) 



Here we use Zubarev's notations.E-tl Within the PM the equations of motion for the GF ( |3.6| ) are projected 
onto the subspace spanned by the operators B ik which leads to the following eigenvalue problem to 



determine and f2 k 



(«). 



r.. _ oMx.. 

l 3 k l J 



) 4 n) W 



where 



-Bi ,k, f 



,4 k }); ({s i>k)J Bt k })=^ 



(3.7) 



(3.8) 



An explicit calculation expresses £jj in terms of the two-site spin-spin correlation functions C r and gives: 



8hirl 



£\2 = £21 = 4r/3 k ^ k , £i 



4fe7T k J k 



£-22 = 



T 
7/2 



(1 - 4 7k ) + | ]T c-^-^C Sl - 



C g (l- 



S2 



+—C S 7 k - 4 - — 



4C g7k + 2C gi+gH dl-51) 



£-23 — £?.2 — 



(3-9) 



where 



- (- cosfc^ + cos k y ) , 7T k 



The projected equation of motion (3.7) defines the mean field Green's functions G$ with three bands 



(i) 
k ' - 



(3) 

il k corresponding to the ZR singlet, the nonbonding oxygen and the triplet bands. 



The mean field spectrum of all three bands is shown in Figs. 1(a) and 1(b) for two different parameter 
sets. In the following we will mainly concentrate on the eigenoperator of the lowest band, the polaron 
eigenoperator B k = i3i ik . The GF of the polaron quasiparticle is defined as: 



G p (k,u) = (B k \Bl 



k/«> • 



(3.10) 



One has to distinguish between the polaron GF defined in ( 3.10 ) and the GF G/((k, ui) which gives the 
number of holes in a unit cell and is introduced as: 



(3.11) 



Please note that according to (2.5) the GF Gh counts not only the holes at oxygen sites, but has a 
contribution at the copper sites also. Below we shall suppose that the intensity given by the photoemission 
experiments can be roughly compared with the hole spectral function which is dictated by G^(k, u>). In 



the mean field approximation G?' (k, u) has the form 



af (k) + a$> (k) 



Gi 0) (k,c) 



E 



- n 



J2 zil) w 



-a 







(3.12) 



which defines the pole strength Z^> (k) of all three bands. If we restrict ourselves to the vicinity of the 

'k 



lowest pole £l k = f2 k ^ the hole GF can be approximated by 
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Gi 0) (k,c)«Z(k)Gf (k,w), (3.13) 

with 

Z (k) = Z« (k) , and = — !— . (3.14) 

F a; — S2k 

Let us first discuss the mean field spectrum of all three bands for the two parameter sets with and 
without direct oxygen-oxygen hopping h (Figs. 1(a) and 1(b)). First of all, one observes a well separated 
singlet band between 1 and 2 eV below the nonbonding one. The singlet-triplet splitting is maximal 
at (7r,7r) and roughly 3 eV (assuming r to be 0.5 eV). That agrees quite reasonably with the splittiM 
between the ZR singlet and triplet states of 3.4 eV that were found in calculations for a CuC>4 cluster £3 
And experimental indications-ibr the triplet state were found in a recent photoemission measurement of 
the substance Ba2Cu304Cl2." It is also interesting to note that the singlet band dispersion obeys nearly 
the antiferromagnetic symmetry (with only small deviations due to the spin rotational invariant ground 
state) but the triplet band not. The nonbonding band receives a dispersion due to h. It has the maximum 
(which would correspond to a binding energy maximum in a spectroscopic measurement) at (n, n). One 
should note, however, that we did not include into our model all the three oxygen 2p orbitals at a given 
site. There exist other nonbonding oxygen bands and especially that one with lowest binding energy at 
(tt, 7r) (seen experimentally in Ref. |37| ) is not present in our model. 

Now, we discuss the lowest spin polaron band fik- Comparison of Figs. 1(a) and 1(b) shows that the 
direct oxygen-oxygen hopping is responsible for the dispersion along the line (0, tt) — (tt, 0). The dispersion 
dependence on frustration in spin subsystem and on temperature was extensively studied in our previous 
workscirEj within an extended basis set. From Fig. 2 we see that the pole strength Z (k) vanishes at the 
r point (with k = (0, 0)). This result is the consequence of a vanishing hybridization between oxygen and 
copper states at this point. Furthermore, the triplet and the nonbonding oxygen bands are degenerate 
there. Consequently, all the orbitals which are incorporated into our model give rise to only one oxygen 
derived peak in the photoemission |Sppctra at the T point. And indeed, the spin polaron is not visible 
in the photoemission spectra there.UEJ Note that the vanishing of the hybridization, and consequently 
the vanishing of the spin polaron is ignored in the course of Emery model reduction to the one-band 
Hamiltonian (see Refs. [l4|,[53] and references therein), and thus the hole spectral weight is nonzero at 
k = (0,0) in the t-J and other jjme-band models. Even the theories for the three-band Emery model 
which start from the Neel ordermLd cannot obtain zero spectral weight at k = (0, 0) since they always 
deal with a linear combination of states with momentum k and k + qo (qo = ("K, tt). 

We note also an apparent similarity of the polaron branch Q±\ < ?fi-fh e present speGfarpp with the calcula- 
tions based on the Neel state and linear spin wave theory for t- JeTBEI and three-bandO~li!3 models. However, 
our results show a slight deviation from the symmetry of the magnetic Brillouin zone ilk 7^ ^k+q due 
to the spin singlet state of the magnetic background. We see from Fig. 1 that the mean field polaron 
bandwidth w ~ r ( w = 2.5r (2.4r) for h = (0.3r) ). That is larger than expected due to the following 
reason: in the limit of small J/t « 1 a state \Bk) = B^Q) for a momentum k near (0,0) or (tt,tt) far 
from the band bottom has an energy Jl max ~ fimin + w ■ This energy may be much greater than that of 
the states 

l^k.q) = ScX-qlS) 

where we defined 

^ = 7^E e ' qR SH- ( 3 - 15 ) 

V R 

The energy of the states |3^k,q) is of the order of Qmin + J f° r k — q near the band bottom. It means that 
|£>k) is unstable with respect to the decay into |3^k,q) states and analogous states, which contain more 
spinwaves. These states contain spin distortions that are situated far from the hole. In order to take them 
into account within the PM we should extend the basis set up to infinity. On the other hand the effect 
of these states can be described by treating the scattering of the spin polaron in terms of irreducible GF, 
also known as the Tserkovnikov technique. As we will see, this repairs the too large polaron bandwidth 
in the simple projection method. 
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4. SCATTERING FOR POLARON GREEN'S FUNCTION 



We treat the polaron operator = B\ k ( |3.5| ) as a candidate for the elementary excitation and calculate 
the corresponding two-time retarded GF G p (k, ui) which is defined in (3.10). The Dyson equation for G p 
has the form 



E(k,w); E(k,w) 



with 



G<,°> = (w - ftk)- 1 ; ft k = [B k ,ff] 
and where we used the irreducible GF 



(4.1) 



(4.2) 



(4.3) 



in accordance with the definition given by Tserkovnikov. The results for G p (k,u>) can be used to 
calculate G/ l (k, w) for low energies 



G h (k,w) MZ(k)G p (k,w) . 



(4.4) 



It must be underlined that Eq. (4.1) coincides only formally with the Dyson equation for the causal Green's 
function. The "self-energy" E(k, ui) in (4.1) has no diagrammatic representation and it is represented 



by more complex Green's functions in the right hand side pf Eq. (4.3). The derivation of Eq. (4.3) and 
its relation to the conventional projection techniqueOE3c3 are represented in Appendix A. As follows 



from (4.1) Re E gives the renormalization of the energy of the polaron B^ and Im E gives its damping 
when Im E(k,o;)| <C |Re E(k, lu)\. In a general case the elementary excitations must be investigated 
self-consistently. The self-energy E(k, uj) is dictated by the interaction of the polaron with the spin 
subsystem, i.e., by the polaron scattering on the spin waves. For this reason the main problem of the 
present technique consists in calculating the irreducible Green's function ( 4.3). 



5. VERTEX CALCULATION 



To determine TZ k we have to calculate the commutator with the Hamiltonian. That gives in detail 
2r/3£ - ) B 1M + 4r/3 k ^kS 2 ,k ^-B 



[5 2 ,k,r] =t 



(31 ~ 3 ' k 



3 Si, k - 2B 2 ,k H 7= 2~1 7k-q G9k-q>l,k,q + 2tk-q>2,k,q) 



£?2,k, J 



£?2.k, h 



4J E 7q\/ |^k-qij,k,q 



= -2hB 2 ,v. 



1 



E 



r l.k,q q r 3,k,q 

Pk-q 



where we introduced the notation 

^,k,q — SqB{ k — q 

In the scattering term of the exchange energy J arises 



(5.1) 



(5.2) 



7 



Y 



J,k,q 



HQ 

— J q 



1 



N 



R 



-i(k-q)R 



(5.3) 



which has a different structure (e Q /3 7 is the antisymmetric tensor). Therefore, it has to be projected 



S2.k, J 



4J- 



1 



E^ 



^k-q^.k.q, 



using 



Finally one finds 



{^■k-q^Vq}) = ~C q ; C q = ^e-^- R ) <S«S£> 

R 



#3,k, H 



4/j-7T k (5k 



I? 



8h% 



l,k 



— i*3,k 



(5.4) 



(5.5) 



(5.6) 



After carrying out the projection (5.4), all terms Rik = j^B^k,-^ 

-Rj,k = Ay-Bj^ic H -= ffij.k.q^j.k^ 



can be represented in unique form 



(5.7) 



with coefficients Ay and 3y,k,q that can be simply derived from (5.1,5.4) and (5.6). One should note that 
the scattering has its origin mainly in the second basis operator i?2.k (3.1) since only coefficients 32j,k,q 
are different from zero. We find from Eq. (5.7): 



where we define 

Rlr = 



= c£\k)R iM = oPCkJAyBAk + R k , 

k = ~m E a ^ ( k )^J,k,q^i,k,q = -J= J2 ( a » (1) ( k ) 9w,k,q^,k-q 



(5.8) 



(5.9) 



We see from Eqs. (4.1,4.3), that the self-energy £(k, ui) accounting for interaction effects is expressed 
through higher or der Green's functions (see Appendix A). One should not ice that the terms linear in 
Sk = £?i jk in Eq. ( [j~8| ) do n ot contribute to the irreducible Green's function (O) for £ (k, uj). The terms 
oc Bj^, ]— 2, 3, i n Eq . ( |5.8| ) are orthogonal to 23k and give nonvanishing contribution to the irreducible 
Green's function (O) only in the energy region of the upper polaron bands. If the lowest polaron band 
may be regarded as isolated (i.e. if it is energetically well separated from other bands), we can neglect 
the interband scattering and retain only the intraband one. This is the case in our problem, since the 
energy gap (1.3t for h = and 1.9r for h — 0.3, see Fig.l) between the two lowest bands is about half 

of the mean field bandwidth of the polaron. Then we project the operator (a\^ (k) gy.k.qi^k-q) onto 

£>k-q and obtain finally 



(5.10) 



where 



T(k,q) = ({af ) (k) fe)k)q B i)k _ q ,ap ) (k-q)S+ k _ q }) = a?> (k) k , q af>* (k - q) . (5.11) 

In the following we will use -Rk instead of !Z k in Eq. (O). Then the lo west order self-energy contribution 
is provided by the first term in the right hand side of the expression (4.3), which is responsible for the 
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scattering of polarons, while the second term leads to higher order corrections (see the discussion in the 
Appendix A). One should men t ion that the operators ^k.q are orthogonal among each other. However, 
unlike the basis operators ( |3.l| , 3.4) only up to terms of order 1/N. Therefore, one may obtain artificial 
terms in the vertex T(k, q) if one carries out the calculation in a changed order than that one presented 

In our previous work§l@ we used only one operator in the PM. The present method generalizes the 
recipe of vertex calculation from Ref. [2^ to the more complex case of several operators in the PM basis 
set. 



6. SELF-CONSISTENT BORN APPROXIMATION 



The self-energy (4.3) with i?k ( 5.1 0| ) provides the opportunity to apply the mode-coupling approxima- 
tion in terms of an independent propagation of the polaron and spin excitations. It consists in the proper 
decoupling procedure for the two time correlation function 



(RuMmlAt')) ~ l^r 2 (k,q)(6 k - q ,,(^_ q!(T (i'))(S- q (t)s q (t')), 



(6.1) 



(s_ q (t)s q (0)„£ eiq ' (r '" 

r,r' 



(S r (t)S r ,(i')> 



(6.2) 



I n th e frameworks of such a mode-coupling approximation the second term in the right hand si de of Eq. 
(4.3) may be neglected (see the discussion in the Appendix A). The first term (Ru]R^)u of (4.2) may 
be expressed in terms of the Fourier transform of the above two-time correlation function ( 6.1). The 
presence of the polaron-polaron correlation function (<Sk- q (£)$k_ q (t')) m ^ ne right hand side of Eq. ( |6.1[ ) 
allows the self-consistent calculation of G p (k, lo). Note that in the case of the t-J model, investigated in 
terms of spinless holes, the analogous decoupling procedure for the irreducible Green's function of the 
form (O) is equivalent to SCBA in a usual diagramm atic technique.0 

The spin-spin correlation function (S_ q (i)S q (t')) in (6.1) is calculated from the spituexcitation Green's 
function D(q, lj) . We treat the spin subsystem in the spherically symmetric approach,t5c3 which is briefly 
described in the Appendix B. The Green's function D(q, w) has the form 



-&J Co 



l"7q 



(6.3) 



col = -32 Jai (Cg/3) (1 - 7q )(2A + 1 + 7q ). 

We neglect the influence of doped holes on the copper spin dynamics and take the spin spectrum pa- 
rameters calculated in Ref. ^5|: A = 0, for T = 0, the vertex correction ot\ = 2.35, the spin excitations 
condensation part m 2 = 0.09. Note that the spin excitation spectrum w q has the magnetic BZ symmetry 
at A = 0, but the Green's function D(q, u>) has the symmetry of the full BZ due to the numerator 

(1 ~ 7q) ' n 

As a result of the decoupling ( |6.1[) we come to the integral equation for the Green's function 

G P (k, u) 1 (6.4) 

u> — ilk — Mk, ^) 

where 

S(k, to) = 1 M '( k ' <tiG P {* - q, w - « q ), (6.5) 
q 
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Af 2 (k,q) = r 2 (k,q) ( 4Cg)(1 7q) . (6.6) 

r(k, q) cpjjpesponds to the bare vertex for the coupling between a spin polaron and a spin wave. It 
is knownEaeSl that this vertex is substantially renormalized for q close to the antiferromagnetic vector 
qo = (it, 7r) . This renormalization is due to the strong interaction of a polaron with the condensation part 
of spin excitations that must be taken into account from the very beginning. As a result, the renormalized 



vertex f (k, q) must be proportional tc 



(q-q ) 2 + I'- 
ll a — > oo in our case of a long range order state of the spin subsystem. Below, this renormalization is 
taken into account empirically by the substitution 



1/2 



L s being the spin-spin correlation length, 



r(k, q ) -» r(k, q ) = r(k, q )^(i + 7q ) . (6.7) 

The introduced vertex correction is proportional to |q — qo| for q close to qo . Let us mention that the 
bare vertex leads to a dramatic decrease of the QP bandwidth. 

7. RESULTS AND DISCUSSION 

In this section we present the results for the low energy part of the one-particle retarded Green's 



function Gh (k, ui) ( 3.11). In our treatment we ignore the scattering of a polaron to the upper bands. 
In this approximation the low energy part of Gh (k, lu) is related with the polaron Green's function 
G p (k, uj) through Eq. ( |4.4| ). We are mostly interested in the spectral function of the polaron or hole GF, 

A p/h (k,cj) = -hm G p/h (k,w) . (7.1) 



respectively, which are given by 



The spectral function of a bare hole Ah is roughly proportional to the intensity in an angle resolved 
photoemission (ARPES) experiment. The self-consistent equation (|6.4|) was solved by means of the 
recursive procedure that provides step by step the coefficients a n ,b n <j*f-ithe continued fraction expansion 
of G p (k, lj). The details of the calculation will be published elsewherecj. The two-dimensional Simpson's 



rule was used for the integration in Eq. (S.5) over 231 points in the irreducible part of the Brillouin zone. 
The prominent feature of the expansion is the fast convergence of coefficients a n , b n to the asymptotic 
behavior which is characterized by a linear dependence on n. The slope for a n is twice as large as the 
slope for b n . This feature allows to use the incomplete gamma function as an analytic terminatoio 
for the continued fraction. We applied the terminator after the calculation of 30 pairs of coefficients. 
Our calculations for various models of strongly correlated electrons indicate that the linear growth of 
continued fraction coefficients seems to be a common feature of such models. 

The solid line in the Figure 1 shows the dispersion e (k) of the lowest peak in the spectral density. 
First, we may see that the overall shape of the dispersion curve has not changed in comparison with the 
mean field result. On the other hand, we may note that the energy renormalization for the top of the 
band is larger than for the bottom resulting in substantial narrowing of the band. Although the band 
width of the polaron in mean field is proportional to the hopping amplitude t, it is reduced to a band 
width proportional to J due to the scattering of the polaron at spin fluctuations which are inherent in 
e (k). Let us note that the polaron mean field energy represents the center of gravity of the polaron 
spectral function A p (k, ui) that is 

/oo 
u)A p (k, oj)dw. (7.2) 
-oo 



On the other hand, we see from the explicit expressions for the Liouvillean matrix elements (|3.8| ) that f^k 
is unambiguously related to the nearest neighbor and next-nearest neighbor static spin-spin cor relation 



functions of the undoped system. The latter is described by the Heisenberg Hamiltonian (2.2) on the 
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square lattice, and t he correlation fvuic .tions are known at present time with a very high accuracy from 
various analyticEScj and numerica£3c3 approaches. That means the asymmetry of f^k (with respect to 
the magnetic Brillouin zone) is not only a consequence of the particular approximations made in the 
present work, but already a property of the polaron spectral function itself. 

Figure 3 gives the polaron and hole spectral densities for three characteristic momenta k to distinguish 
three regimes. In addition, we show the real and imaginary part s of the self-energy £ (k, ui) which are 



approximately identical for polaron or hole GF according to Eq. (4.4). Let us recall that we have made 
the change k — > k' = k — qo in the final results, so that the quasimomentajja the figures correspond to the 
actual Brillouin zone of the CuC>2 plane, in contrast to our previous work.ES We find a sharp quasiparticle 
peak at the bottom of the spectra for k a = (ir/2,ir/2) and kb = (77r/20, 77r/20). The position of the 
quasiparticle peak corresponds to the condition Re G~ 1 (k, ui) = 0, i.e. the point where we have a crossing 
of the functions u> — J7k and Re £ (k, ui), see Figs. 3(a) and 3(b). Quite a different shape has the spectral 
density A p (k c , uj) at k c = (0,0) . Its lowest peak is not related to the zero of Re G" 1 (k, uj) and is due 
to abrupt changes in the self-energy. That leads to a large imaginary part of £ (k, uj) at the position of 
the lowest peak in the spectral function. We may say that there is no quasiparticle excitation for this 
quasimomentum value. 

Figure 4 shows the difference between the quasiparticle excitations at points k a and kf, more in detail. 
The first one corresponds to a quasiparticle with infinite lifetime, i.e. to an isolated pole of the Green's 
function at the real axis of the complex energy plane. For k(, the situation is different, a very small, 
but finite, imaginary part of the self-energy is present in the region of the peak. The shape of the peak 
deviates from Lorentzian. Therefore, such an excitation can be regarded as a quasiparticle with finite 
lifetime. 

Table 1 shows the position of the lowest peak, e (k) , the area under the peak for polaron W p (k) and 



for hole spectral function, Wh (k) = Z (k) W p (k) (see (4.4)), and the value of the imaginary part of the 
self-energy, — Im £ [k, e (k)] for k values which vary along the diagonal of the Brillouin zone. The rows 
which are marked by a star do not represent a quasiparticle excitation like at k c and the filled or open 
circles correspond to momenta with quasiparticles of infinite or finite lifetime, respectively. 

So, we may distinguish three qualitatively different kinds of spectral function A p (k, uj) behavior: 

(i) The lowest peak tends to a delta function W p 5 [uj — e (k)], when Im ui — > (see Fig.4(a)). We have 
a quasiparticle, characterized by the excitation energy e (k) = + Re £ [k, e (k)] with infinite lifetime 
T\ , (1/tj) = — Im £ [k, £ (k)] = 0. These k-points are situated near the band minimum k m j n . They are 
marked by solid circles at Fig. 1 and in Table 1. 

(ii) The lowest peak has approximate Lorentzian form (see Fig. 4(b)) 

A p (k, u) = --Im G p (k, uj) » W p - + A mcoh (7.3) 

7T [w-e(k)] +{l/n) 

and corresponds to a quasiparticle with energy e (k) = + Re £ [k, e (k)], and finite lifetime (I/17) sa 
— Im £ [k, e (k)] <C |e (k) — e (k m i n )|. The corresponding k- values are marked by open circles at Fig. 1 
and in Table 1. 

(iii) The spectral density is completely incoherent and is of the same order of magnitude as the imagi- 
nary part of the self-energy at all uj (see Fig. 3(c)). These k-points are ma rked by stars. 



In order to understand the nature of the polaron damping (1/t; in (7.3)), let us consider a particular 
k- value and est imate the imaginary part of the self-energy at a particular uj. Changing the summation 
variable in (|3.5| ) we can calculate the imaginary part of the self-energy as follows 

-Im£(k, W ) =7rl^M 2 (k,k-q'M P (q>-^ k _ q . (7.4) 
q' 

Now we suppose that the lowest peak at all q' = k — q may be described approximately by Lorentzian 



( |7.3[ ). Then we see from Eq. (7.4) that large damping, characterized by a large value of — Im £(k, ui) may 
only arise when the condition 

uj - e (q') - lu^ = (7.5) 
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is satisfied for some values of q' (see also Ref. 14). On the contrary, if Eq. (7.5) does not hold even at 
the lowest edge of the spectral function ui — e (k), the damping is absent. In Fig. 5 we plot e (q') together 
with the curve e (k) — CL>k_ q ' for various values of k (for h = 0.3r and where we choose only one direction 
for both k, and q' ). We may see that for k m ; n =(7r/2, tt/2) the condition ( |7.5| ) holds only for the trivial 
values q = k — q' = (0,0) and q = q , where, however, the vertex M 2 (k, q) vanishes. It demonstrates 
that it is impossible to scatter from k min into other states. The same situation occurs for some finite 
region of k values around the band bottom. These are exactly those momenta with a quasiparticle of 
infinite lifetime. Different is the situation for k=(0,0) and k =(47r/20, 47r/20) where several nontrivial 
intersections take place. 

Figure 6 shows the low energy part of the hole spectral density Ah (k, lu) for k along high symmetry 
directions. Let us recall that the intensity of the lowest peak o f Ah (k, uj) is governed by two circumstances: 
the k-dependence of the mean field residue Z (k) (see Eq. ( 3.12 ) and the discussion about Fig. 2) and 
the intensity of the lowest polaron peak of A p (k, oj). One can clearly observe a well defined quasiparticle 
peak near the bottom of the spectrum at (ir/2, n/2), everywhere along the line (0, 0) and also near 

to (-7T, 0). At the T-point (0, 0), there is no hole spectral density due to the vanishing residue for the hole 
Green's function in the lowest mean field polaron band. At the same time one observes intensity at (n, tt). 
A clear asymmetry of the peak intensity is seen along the diagonal of the BZ. The abrupt drop of the 
peak intensity in the region (7r/2,7r/2) -(7r, tt) is related to the strong polaron damping there. In other 
words, most of the mean field hole spectral weight {Z{k) in Fig. 2) goes into the incoherent part of the 
spectrum (situated at higher energies) due to the strong coupling of the polaron with spin excitations. 
An analogous behavior of the quasiparticle peak intensity was obtained in Ref. |l5| within a variational 
ansatz for the t-J model. 

Let us compare the results with those of the usual SCBA of spinless holes in the pure t-J model. 
There, one finds quasiparticles with infinite lifetime and finite weight everywhere in k-space. Addition- 
ally, the spectral function has the symmetry of the magnetic BZ. Introducing additional hopping terms 
(corresponding to direct oxygen-oxygen hopping in the Emery model) leads to a scattering mechanism 
such that the upper parts of the spectrum loose their quasiparticle character.E3o On the contrary, in 
our approach, the damping of the spin polaron is already present without direct oxygen-oxygen hopping. 
If we consider only the dispersion relation e(k) we observe a remarkable similarity between the present 
calculation and earlier results using the SCBA of spinless holes or other methods. In fact, the devia- 
tions from the magnetic BZ symmetry in Fig. 1 are not very large. On the other hand, we find some 
qualitatively new features of the spectral function which are absent in the numerical results of Refs. ||,(7| 
Besides the strong polaron scattering away from the band bottom we note especially: (i) the absence of 
the polaron quasiparticle at the T -point and (ii) the asymmetry of the peak intensity with respect to 
the magnetic BZ. 

Our results (Fig. 6) can also be summarized in such a way that a well pronounced polaron quasiparticle 
peak exists only around the bottom of the band. In the diction (tt/2, tt/2)-(tt, 0) it is more clearly seen 
than perpendicular to it. The recent experimental findingEJ shows, indeed, that the Zhang-Rice singlet 
can be observed in Sr2Cu02Cl 2 only in a similar region of the BZ. And also the higher peak intensity 
going from (7r/2,7r/2) to the T -point in comparison with the opposite direction is in agreemejit jyith 
experiment. So we see that those details in which our calculation differs from the standard oneaoETLl are 
essential for a better understanding of the experiment. 

Let us emphasize that for a detailed comparison to a real experiment there are a considerable number of 
complications which have not been taken into account in the present work: First, in general a photoholc 
can be generated also on copper, and there will be (in general k-dependent) interference between the 
phot oholcs created on different atoms. An operator cj^ +a (|2.5| ) that enters the spin-fermion Hamiltonian 
( |2.l| ) creates a hole in a state that is an antibonding (in hole notation) combination of oxygen p a state 
and singly occupied copper d x 2_ y 2 state at adjacent sites. In a more realistic theory one should deal with 
a complete set of bare oxygen and copper hole creation operators. This may be expected to have some 
influence on the spectral weight. Next, one must take into account interference terms like < p x \Py >, 
because the photoemission operator is in general a linear combination of p x and p y . Thereby the relative 
phase between p x and p y depends on the momentum transfer, and is in general even different whether 
one is in the 1st or 2nd Brillouin zone etc. Clearly this will give additional momentum dependence of 
the weight. Moreover, it is not even a priori clear that p x and p y have equal weight in the PES-operator. 
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Rather this depends in a relatively complicated way, e.g., on the polarization of the incoming X-ray 
photons. It is in fact well known experimentally that the spectra depend-auite sensitively on the X- 
ray polarization in several oxychlorides like S^CuC^C^ or Ba2Cu3O4Cl2Ke.il At the present stage of 
many-body theory it is impossible to take all these complications into account. Thus, our comparison 
with experiment has only qualitative character. Nevertheless, we point out that according to our results, 
many features of the experimental spectra may originate from the peculiarities of the correlated ground 
state of the material rather than from the above complications. 

We may consider now what happens if we begin to dope our system. First, we suppose that for 
extremely small doping the Fermi surface will consists of hole pockets centered at (tt/2, ir/2). The small 
spectral weight of the quasiparticle poles in this region means that the surface enclosed inside the Fermi 
surface may be much larger than the number of holes. Note that due to the strong weight asymmetry the 
experimental observation of the Fermi surface parts facing M point with k = (71", tt) point may be very 
hard to do. In Ref. [l5]the observation of weak features in photoemission spectra in this region of k-space 
was reported, but the authors can not unambiguously interpret them as Fermi surface cuts. With the 
increase of doped holes the antiferromaenjetic correlations in the spin subsystem are weakened resulting 
in a deformation of the hole dispersion. E£lo The minimum of the hole spectrum should shift to the M 
point and a 'large' Fermi surface develops. An analogous scenario was also developed in Ref. [lO]. 

8. CONCLUSION 

We proposed a new method to calculate the spectral function of the Emery model in its spin-fermion 
form combining the Mori-Zwanzig projection method with a representation of the self-energy by irre- 
ducible GF. Self-energy effects reduce the mean field bandwidth of the spin polaron to a realistic band- 
width of the order of 2J, but do not change the overall shape of the dispersion. Direct oxygen-oxygen 
hopping leads to a more isotropic band minimum around (7r/2,7r/2). The quasiparticle weight is maxi- 
mal (about 0.25) around the minimum of the dispersion. Despite some similarities- wie observed several 
new features of the spectral function which are not present in the usual SCBAaLffil for spinless holes 
coupled to spinons in the pure t-J model, namely the absence of quasiparticle weight at the T point, 
the strong damping away from the minimum and the asymmetry with respect to th e . an tiferromagnetic 
BZ (where the asymmetry is present, however, in other approaches to the t-J model)t3ll3. Qualitatively, 
these features can be observed in the ARPES experiments on Sr2Cu02Cl2- We studied in detail possible 
damping processes of the spin polaron and distinguished three regions in k-space where the spin polaron 
exists with infinite lifetime, finite lifetime or does not exist at all. 
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Appendix A 



Here we give the derivation of the Tserkovnikov expression for the self-energy (4.3) and compare it 



with the conventional projection technique expression. We use the equation of motion for the retarded 
Green's function G p (k, u) = (s k | B£) Q 



co(B k \Bl 



n k | bA 

I u 



1 + (Bk I ft£ 



(A.l) 



Multiplying the equation for the higher order GF 

u (n k \ St)^ = ({ll k ,Bi}) + ( R k I R k 

by ^B k | B k ) , and noting that ^|7£ k , B^ = £l k , we have 



wfaklBl) (B k \Bt) = ft k (B k |B, 



In the left-hand side we use the second equality (A.l) to obtain 



K*\Bt 



n k (B k \Bl 



n k \nt) (B k \st 



n k I b. 



B k I ^\ 



and substitute it into first equality (AA). Finally Eq. ( |A.l|) takes the form 

[w - fi k - S (k, w)] Gp (k, w) = 1 



(A.2) 



(A.3) 



(A.4) 



with the self energy given by (4.3) 



In order to clarify the physical meaning of this expression let us compare it with the close expression for 
Spt (k, to) which follows from the projection technique (see Refs. |3l] , |3^ , |38| ) . Taking the usual notations 
for the scalar product in the operator space, the projection operators P and Q and the Liouvillean 
superoperator C 



(A\B) = ({A,B}) , P =| 4)(B k |, Q = l-P, CA = [H, A] , AC = [A,H] 
Green's function and the self-energy may be written as 



(A.5) 



G p (k,w) = (B k 



— — BA 



1 



si). 



E PT (k,w) = (B k | CQ 



1 



to — CQ 



CBl) = (B k CQ 



uj-C 1 ~ k 

| 1 

to — QCQ 



QCBt), 



(A.6) 



The Tserkovnikov expression (4.3) has the following form in the projection technique notation 

1 1 „+, 1 ,„ , 1 



E T (k,w) = (B k £ 



(B k £ 



(B k £Q | — l —QCB\.) - (B k CQ 

LU — L, 



rB 



£ (Bk I zj^BD 



(Bk I ^B k ) 
^— K(Bk I -^Q^t) 



(A.7) 



In the last equality we excluded the terms linear in B k . Now we see, that the first term in Ep (k, to) 
describes the propagation of a "fluctuation" B k CQ in the full Liouvillean space, and the second counter- 
term eliminates items corresponding to the propagation in the subspace spanned by B k , i.e., it makes 
the same job as the projectors Q surrounding C in the denominator of Epp (k, to) (A. 6). Using the 
known rules of matrix algebra it is in fact easy to see that Epp (k, to) ( |A.6| ) is equivalent to Ep (k, to) 
(A.7) . W hen we consider only intraband scattering we approximate B k CQ ~ i? k , with 7? k , given by 
Eq. (|5.10|) . Then outset of approximations is equivalent to the mode-coupling approximation in the 
projection technique^ 
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Appendix B 

Now we review briefly the Kondo-Yamaji theoryi! for an isotropic Heisenberg spin-half antiferromagnet 
on a square latticeO'Ej 

In the absence of holes our Hamiltonian reduces to the spin Hamiltonian J (^^). We search the 
retarded spin-spin Green's function -D (q, uS) = (5^ | S-q)^ taking into account that (S^) = and 



oa CP 
^R^R+r 



i fia0(j r m t ne adopted spherically symmetric approach. Then we have the following 



equations of motion for the spin Green's functions 



to 



n.m 



s 



(B.l) 



zJ\C g \ (1 - 7q ) + -z J 2 (1 - 7q ) w (S£ I S« q )^ + 



1 Pt~& 

t2 \ ' iq(n— m) \ * / q/3 q/3 no , rr/3 c/3 m_ ( q/3 qP i c/ 5 c/3wa i ™ 

J / y c / y y J n J n+g"n+g- gi t J n+g"n+g 1 "n l"n"n+g 1 ' J n+g-gi J n/ J n+g I °m 



E 

2j|C g |(l- 7q 



jV-ICglarCl-T,) (2A + 1 + 7,) 



(B.2) 



2 = 4 is the coordination number. The last approximate equality was obtained by the unique decoupling 
procedure for three-site (on different sites) spin-operators: 



(B.3) 



Here the parameters ai or ui are used for the first and the second nearest neighbors n and m accordingly. 
They are regarded as vertex corrections which have to be introduced in order to prevent the violation of 
the sum rule of the correlation function. The energy gap parameter equals 



2A 



"2 J2 g , Sl c g+gi a 2 -l 1 



z- 1 



iiz 2 ICg 



ai 2z\Ca 



(B.4) 



The substitution of Eq. ( |B.2j ) into Eq. ( |B.l| ) provides D (q, uj) in the form ( p.3| ) . The static correlation 
functions C r = (SrSr +i .) may be expressed through D (q, uj): 



N ^ 



-?qr 



^ MmJ}( y )/ffl 4^e- 



exp(uVT)- 1 



AT 



~zJCg(l-7 q ) 



coth (w q /2T) 



(B.5) 



For r = 0, Co = 3/4, and the e quat ion (B.5) gives the sum rule for D (q, a;). Thus we obtain the set 
of self-consistent equations ( |6.3|) ,( [B~5|) . Detailed considerations show that it is overdeterminedo and one 
extra condition may be imposed to solve it. This condition is the relation between a\ and ct2 given below. 
In a macroscopic system (N — * oo) long-range-order is present for T = 0. We have A = 0, 



uj„ = zJ (2ai ICgl /3) 1/2 1/1-72 m 2 = lim c 4qor C r ^ 0. 
The set (B.5) reduces to 



-y 

N ^ 



q#q 



ZJC S (1 -7 q ) 



e- iqor m 2 . 



(B.6) 
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Shimahara and TakadaEa explored various kinds of extra conditions to determine the parameter's values 
and found m to be the most sensitive quantity. They proposed to fix r a — (a% — 1) / (a^ — 1) = 0.82579 
for all temperatures which gives m = 0.3 for T = 0. Then the values of ground state energy and uniform 
spin susceptibility agree with other theories and Monte-Carlo simulations within few percents. For finite 
temperature the KY theory leads to the absence of LRO and provides a spin susceptibility which has a 
satisfactory behavior over the whole temperature region: it has a zero temperature derivative at T = 0, 
it exhibits a peak around T ~ J and reproduces the high temperature expansion for T > 1.6J. In the 
absence of LRO, the low lying excitations are essentially spin waves propagating in a short-range order 
with a temperature dependent correlation length. On the other hand, if there exists LRO, the usual spin 
wave excitations are recovered. 
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FIG. 1. The dispersion of the lowest QP band e (k) and the mean field dispersions Q,^ along high symmetry 
lines in the Brillouin zone for J = 0.2 and r = 1 without (a, h = 0) and with (b, h — 0.3) direct oxygen-oxygen 
hopping. Filled and open circles mark quasiparticles with infinite and finite lifetime, respectively. A star cor- 
responds to a lowest peak in the spectral density which is strongly overdamped. The points in k space mean: 
r = (0,0), X = (tt,0), M = (7r,7r) and X' = (0,tt). 



FIG. 2. The QP weights Z^' (k) of all three bands in the projection method (PM) and the area under the 
lowest peak of the hole spectral function Wh (k) calculating the irreducible GF of the self-energy in self-consistent 
Born approximation (SCBA) along high symmetry lines in the Brillouin zone for J = 0.2 and r = 1 without (a, 
h — 0) and with (b, h = 0.3) direct oxygen-oxygen hopping. 



FIG. 3. Polaron and hole spectral densities A p / h (k, uj) around the lowest peak for the momenta: a) 
k = (7r/2, 7t/2), b) k = 7(7r/20, 7r/20) and c) k = (0, 0). Also shown are the imaginary part of the self-energy and 
the crossing of the curves uj — il^ and Re E (k, a>). The spectral functions are broadened by a small imaginary 
part n — 0.005 in uj. 



FIG. 4. Polaron spectral density and self-energy for: a) k = (tt/2, vt/2) and b) k = 7(7r/20, 7r/20) to distinguish 
a quasiparticle with infinite lifetime (a) and a finite one (b). 



FIG. 5. Crossing of the curves e(q') and e(k) — <^k- q ' f° r three different values of k along the diagonal of the 
BZ for q' = (q' x ,q x ): there are only trivial crossing points (q = k — q' equal to T or M) for k = (n/2, 7r/2) but 
nontrivial crossings for k = (4-7r/20, 4-7r/20) and k = (0, 0). 



FIG. 6. Contour plot of hole spectral function Ah (k, u>) for for J — 0.2, r = 1, h = 0.3 and n = 0.01 along 
high symmetry lines in the BZ a) (0, 0)-(tt, tt), b) (tt, 0)-(0, tt) and c) (0, 0)-(tt, 0). 
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TABLE I. Position of the lowest energy peak e(k), the areas under the peak W p (k) and Wh(k) and the 
imaginary part of the self-energy for J/r = 0.2, h/r = 0.3 and k = n(7r/20, 7r/20) 
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